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We study the static structure factor of the fractional Chern insulator state. The averaged over 
Brillouin zone Fubini-Study (quantum distance) metric of the underlying non-interacting quantum 
system enters the quadratic form with small momenta in the expansion of the single particle part of 
the static structure factor. The form corresponds to the plasmonic part in the Laughlin case, and 
thus we find that the averaged over Brillouin zone Fubini-Study metric plays the role of "Landau 
level" metric in the framework of the geometrical description of fractional quantum Hall systems 
[F.D.M. Haldane, PRL 107, 116801 (2011)] and an effective (continuum) description of fractional 
Chern insulators. Assuming Laughlin-like correlations in the two body part of the static structure 
factor, we discuss the properties of the static structure factor in the long-distance limit, and an- 
alyze the conditions that have to be satisfied in order for the FQHE scenario [S.M. Girvin, A.H. 
MacDonald, and P.M. Platzman, Phys. Rev. B 32, 8458 (1985)] to occur. We discuss the relation- 
ship of the (averaged) Fubini-Study metric and Berry curvature, and illustrate their influence and 
correspondence in Haldane model based fractional Chern insulators. 



I. INTRODUCTION 

Chern insulators (CIs) ilj exhibit integer quantum Hall 
effect (IQHE) conductance quantization in the absence 
of the magnetic field due to non-trivial filled band struc- 
ture with non-zero topological Chern number. Fractional 
Chern insulator (FCI) is the name for CI with a par- 
tially filled band (which is akin to a Landau level) in the 
presence of strong interactions, which exhibit fractional 
quantum Hall effect (FQHE) conductance quantization. 

What might be defined as an ultimate goal in the 
context of the FCI physics would be the understand- 
ing of the mechanism of creation of FCI states in order 
to be able to suggest the most convenient experimen- 
tal settings, whether in cold atoms, crystal (solid-state) 
physics, or graphene structures for their realization. In- 
teracting topological insulator |SH3 physics, which FCI 
is a time-reversal-symmetry-broken representative, is one 
of the major topics of the current research in the field 
of strongly correlated systems. FCI have been reported 
in many models both for fermionic systems |1H11 IH] and 
bosonic cold atomic models [HI HO] ■ Still there is no com- 
plete understanding why some model crystal systems are 
more convenient than others for particular FCI states. 
One way to resolve this problem would be to study the 
underlying quantum geometry of the crystal system and 
its influence on the gap function of FCI in the scope of the 
usual approximation in the FQHE physics: single mode 
approximation (SMA) [TT] . 

In practice finding favorable conditions for FCI means 
achieving the understanding of the interplay of the statis- 
tics of the underlying particles, interactions, and the 
background lattice that produces the band with non- 
trivial topology characterized by a non-zero Chern num- 
ber. We will focus our attention on the influence of the 



background lattice, and the structure of its particular 
band in which the FCI physics takes place. The band 
structure is characterized by two tensors that provide the 
description of the evolution of the quantum - mechanical 
state as it changes with the change of the Bloch mo- 
mentum of the lattice: the Fubini-Study (FS) metric (or 
quantum distance) and the Berry curvature. The ques- 
tion of the role of the FS metric in the context of FCI was 
first raised in Ref. . A recent work[T3] showed how the 
FS metric is related and can be measured via the current 
noise spectrum. While there is a general understanding 
that an almost fiat Berry curvature favors the emergence 
of a FCI phase (mimicking the constant magnetic field 
of FQHE), there is a little or no understanding what is 
the infiuence of the FS metric. We know that there is no 
quantization of its value when integrated over the Bloch 
momentum phase space - BrilUouin zone, like in the case 
of Berry curvature (that produces the Chern number). 
Thus quite generally its averaged value is dependent on 
the lattice parameters, and we cannot reduce the situa- 
tion to the one of a single Landau level (FQHE) in which 
the relationship between, what we may identify as. Berry 
curvature and FS metric is fixed and not model depen- 
dent. 

On the other hand we want to point out that there 
is the largely unexplored problem and question on the 
status of SMA for FCI. In the search for favorable condi- 
tions for FCI this is the most likely tool that is available if 
we want to understand the parameter (FS metric) depen- 
dence of the measure of the stability of the FCI state - its 
gap function. In this work we will pay special attention 
to the static structure factor (SSF), the norm of the ap- 
proximate state, although "f(k)" (the oscillator strength) 
function, i.e. the expectation value of the energy of the 
system with respect to this state is also important and 
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needed. In the literature we find a general statement 
on the nature of the expansion of the f(k) function for 
FCI; the leading term in the long distance expansion is 
quadratic in the low momentum due to the fluctuating 
Berry curvature in FCI. This necessarily means, in the 
scope of the SMA, that the SSF has the same leading 
behavior. We want to point out that the derivation in 
the same reference (14j was done with the usual averag- 
ing and long-wavelength limit procedures [TStiTT] in FCI, 
but which were not assumed or applied in the original 
work in Ref. [Tl] on the SMA in the context of FQHE. 
The behavior that was described in Ref. [TT] is the lead- 
ing behavior of the SSF and f(k) as quartic in small mo- 
menta, and this is, we may say, a hallmark of the FQHE 
behavior. 

Thus it is natural to ask the following question: 
Whether, in a mean field picture at least we can talk 
about the scenario of Ref. [TT] for FCI, or lower order 
terms in the SSF and f(k) that depend on the fluctu- 
ations in the Berry curvature and FS metric will de- 
termine the gap function. In this work we will study 
the SSF for FCI Laughlin-like states. In doing this one 
may resort to a numerical calculation of the SSF of FCI, 
which undoubtedly will shed most light on the nature 
of the expansion of the SSF - whether or not we have 
the quadratic term. Due to the approximations made 
and complicated expressions in Ref. [T^ that needs to 
be checked for concrete FCI states. On the other hand, 
as will be described in this work, one may build a mean 
field picture (valid at least in the constant curvature and 
constant metric case) based on analytical considerations 
in the long-distance limit. By analyzing the single par- 
ticle physics of the band with non-trivial Chern number 
and its consequences on the SSF of FCI, we can find out 
that the effective ( "plasmonic" ) density-density force in 
the long-distance limit is inversely proportional to the 
averaged over Brillouin zone (BZ) (non-zero) diagonal 
element (s) of the FS metric. Thus the averaged over 
BZ FS metric plays the role of "Landau level" metric in 
the framework of the geometrical description of fractional 
quantum Hall systems |18| and an effective description of 
FCIs. In the effective description of the non-interacting 
CI band the density-density force will ensure that the 
long-distance fluctuations are suppressed - that the non- 
interacting system is gapped, and the density is uniform 
in the long-distance limit. By applying the Feenberg for- 
mula |19fl21j interacting quantum liquid we can track 
down the influence of the FS metric on the ensuing coef- 
ficient of the quartic term of the SSF which will generalize 
the expression of Ref. [H] to the FCI case. 

Thus the problem of calculating the SSF in the context 
of FCIs possesses many difficulties, and makes transpar- 
ent the very nature of the FCIs' complex background 
i.e. fluctuating Berry curvature and fluctuating FS met- 
ric. The geometric picture of FQHE systems proposed in 
Ref. [E] , in the context of FCIs may find the most non- 
trivial realization where also extrinsic part of the metric 
degree of freedom (not the one related to interactions) is 



non-trivial. In this work, a definite answer for the SSF 
behavior in the long-distance limit will not be given, but 
we will analyze the most relevant properties in the same 
limit by analytical means. It will be explained which 
conditions have to be fulfilled in order for the FQHE sce- 
nario (Ref. [11 ) to occur in these systems. To illustrate 
our main conclusions we analyzed the underlying quan- 
tum geometry of the Haldane model p[j and the pres- 
ence of bosonic FCI states in the phase diagram of the 
same model with on-site interaction (repulsion) only. We 
choose bosons in order to minimize the influence of the 
symmetry of the underlying lattice. The paper is orga- 
nized as follows. In Section |ll] the derivation of SSF in 
the context of FQHE is reviewed. In Section [TTT| the role 
of the FS metric in the context of FCIs and FCI SSF is 
identified. In the scope of a mean field approximation to- 
gether with the assumption of the Laughlin-Jastrow cor- 
relations, the long-distance limit of the SSF for Abelian 
FCI states is calculated in Section |lVj The possibility 
for the same scenario of Ref. [TT] in the context of FCI 
states is discussed in the same section. In Section |V| 
the background degrees of freedom, the FS metric and 
Berry curvature, for the Haldane model based FCIs are 
analyzed. In the same section, Section [V| the phase dia- 
gram of the interacting system of bosons (that live on the 
lattice defined by the Haldane model) is presented, and 
the identification of possible FCI regions and a compari- 
son with the calculated background properties was made 
possible. 

II. THE STATIC STRUCTURE FACTOR AND 
LAUGHLIN CASE 

A. Static structure factor 

We deflne static structure factor (SSF) as 

s{q) = ^{^\p-,P,\^), (1) 

where V is the volume of the system, 1 4') is a normalized 
many-body wave function, and 

y dr exp{-iq • r}p(r), (2) 

where 

JV 

p(r)=^<5(2)(r-r,0, (3) 

i=l 

is the density operator of the system of N particles with 
coordinates {r^; i = 1, . . . , N}. If we introduce the radial 
distribution function, 

g(|ri-r2|) ^ J d\s--- J dr^|*(ri,...,rA.)|2, 

(4) 
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where n = N /V is the averaged density, we can rewrite 
Eq.[l]as 



s{q) = n + J dr g{r) exp{— iqr} 



(5) 



We wiU cah the first term in s{q) single particle part and 
the second one two-body (correlation) part. 

In the second quantization language with creation and 
annihilation operators in the coordinate space, ^^(r) and 
#(r), 



where so(fc) represents the contribution from "irreducible 
diagrams". In the lowest order in m, so(q) ~ n = l/27rTO. 
We set the magnetic length to one, Ib = I- The single 
particle contribution, n is contained in Eq. [9j 

In the lowest Landau level (LLL) we can define the 
projected density. 



Pi 



d 



d"z*T(^)exp{ig^}exp{i^}*(z), (10) 



where ^'(z) are second quantized operators in the LLL, 



[*(r),^^(r')]± =(5(2)(r 



(6) 



where + sign denotes a commutator in the case of bosons, 
and — sign denotes an anticomniutator in the case of 
fermions, we can easily work out the corresponding ex- 
pression for s{q): 



5('Z) = ^(*I / dr 4't(r)^(r) |vl/) 



+ — I dr I dr' (*|^^(r)'i'^(r')*(r')^'(r)|*) x 



exp{— iq(r — r')}. 



(7) 



In other words, by using the basic algebra of particle 
operators i.e. (anti) commuting relation in Eq. 6] we can 
recover separate single part (first term of Eq. ^ ) and the 
two-body part (second term of Eq. [7|. 



B. Quantum Hall case - Laughlin case 

We would like to calculate the same quantity SSF, de- 
fined in Eq.[l] where j^I') is the normalized Laughlin wave 
function. With respect to Ref. |^ , we use the definition 
of SSF given by Eq. [l] in which we divide the density- 
density correlator by volume (V) instead of the number 
of particles (N). The definition used in Ref. [55] corre- 
sponds to the norm of the model state used in SMA. In 
the following we will summarize the results of weak cou- 
pling plasma approach as described in Refs. [5T] and [55] . 

If we assume weak coupling in the Laughlin plasma 
approach [24] 125] we can study expansions of expectation 
values as in Eq. [T] in terms of the 2D Coulomb plasma 
interaction, 



v(q) = — n2", 



(8) 



where positive integer m is connected to the filling factor 
v of the quantum Hall system as ly = 1/m. Because of 
the expected screening of the 2D plasma the contribu- 
tions in the weak coupling perturbative approach can be 
organized as in 



1 



ai 



1=0 



'exp{--|2p}, 



(11) 



z ~ X + iy is the complex 2D coordinate and [d;, dj„]± — 
Si^rn- The derivatives in Eq. 10 act only on the holomor- 



phic part (dependent only on z) of ^'(z). We have 



P-qPq 



d2z'[exp{-ig*^}exp{-i^}vl/t(/)]vl/(z') 



d 



x¥{z) exp{ig^}exp{i^}*(z), 
oz 2 



(12) 



Similarly to the unprojected case, Eqs. [6] and [7] we can 
use the (anti)commutation relations of projected single 
particle operators, ^"(2), to get the division in the single 
and two-body part of SSF. The single particle part of 
the SSF corresponds to the diagonal contribution either 
in momentum or coordinate space that refers to a single 
particle. We wiU denote by {p-qPq)\ single J the single part 
of the projected SSF. Due to 



[vl/(z'),^n^)]± = ^exp{^}exp{-M!}exp{-^}, 



(13) 



i.e. equality of the (anti)commutator to the LLL delta 
fmiction, the single particle part of the projected SSF is 



(P-qPq) I single = ^ exp{ — }. 



(14) 



The two-particle correlations (i.e. when z' ^ z) stay 
the same (examine Eqs. 12 and [t] with Laughlin wave 
function in that case) and therefore the expression for 
the projected SSF is 



s(g) = nexp{ —} + — 



so(g) 



(15) 



The important correction to the sq {q) ~ n approximation 
is made of a bubble diagram in the next order in m; the 
interaction lines (two of them) are screened. Therefore 
the contribution is 



saiq) 



l-v{q)so{q)' 



(9) 



1 



d^fc 



V;ff(q-k)T4ff(k), 



(16) 
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where 1/2 is a symmetry factor and 



Viq) 



1 - nV^(q) 
The correction for q = is 

m 

do = n— 



where Up = (7^7p) is the occupation of the Bloch mo- 
mentum p in the many-body FCI state. 

(17) After a few steps, which are described in Appendix A, 
the expansion in small momentum q of the single part of 
the unprojected SSF for the FCI state is 

(18) {P-q Pq)\singic = n- q^qJng^f + -^^ndkgfj" + o(q ) 

(25) 



We checked that the |qp correction to 6q in so{q) ~ n+Sq where we assumed summations over repeated indices, 
will not change the resulting behavior due to the inclusion 



of SQ{q) ~ n + 6o in Eq. 15 



(m-l) 4 
3 111 



(19) 



If we assume that the analytical continuation for m > 1 
is valid we recovered the formula of Ref . [22] that follows 
from the 2D plasma compressibility rule. 



III. FCI CASE - SINGLE PARTICLE 
CONTRIBUTION 

We start with a lattice system in which Bloch Hamil- 
tonian H{p) is diagonal with eigenvalues, emip), where 
m = l,...,r denotes the band index, and eigenstates, 
where a = 1, . . . ,r denotes the orbital index in a 
unite cell, and p is the Bloch momentum, i.e. 



(20) 



In order to study the partially filled band with Chern 
number equal to one in which a FCI state occurs, we 
confine our description to that band, and drop the in- 
dex m in the following. In the case of FCI we take the 
projected density to a single band to be defined as in 
Ref. HZ] i.e. 



E* t 

P 



(21) 



where the summation on the repeated orbital Greek in- 
dex (a) is assumed, p and q are Bloch momenta, and jp 
are normal mode operators. 
Therefore 

P-qPq = 

^ ] ^a,pi'^a,pi-qU^p^Ua^P2+q X 



and because 



we have 



[^pi-giTpali ~ ^Pi-q,l 



(P—ijPg) I single 



(22) 
(23) 

(24) 



9i'Uai,pU^^p'Uc,2.pdjU^^p — djUaj^,pU^_^pUa2.pdiUg^^p] 

(26) 

is the Fubini-Study metric, see, for example, Ref. [T^ for 
the definition, and 



■p 



and 



i:,d,gff{p) 



(27) 



(28) 



are averages over the whole BZ. If we assume Up — 
constant, which we expect to hold at least approximately 



in the FCI state, see Section IV B (Eq. |48|) for an expla- 



nation of this point, we can substitute the averages over 
occupation number to the ones over Brillouin zone. Be- 
cause of the periodicity in the A;-space we have that the 
averaged derivatives over the metric (Eq. 28) are equal 
to zero in the case of FCIs. Also we can argue that, as- 
suming the inversion symmetry, the expression in Eq. [28] 
is zero. Thus 



{P-qPq)\ 



single 



n - qtqjugff 



o{q' 



(29) 



Comparing with the Laughlin case, Eqs. [14] and [T5| we 
expect that the quadratic term corresponds to the plas- 
monic part that has to be canceled in the projected (to 
the band) SSF if the FQHE scenario occurs in the context 
of FCI. 

Continuing the analysis of the single part to the quartic 
order we have 



-,FS 



(p-qPq) Isingic = n- qiqjU gl-" + qiqjqkqin n^ju 

+ higher order terms (30) 



where 



-^ [ua,pdkdidjdiu*^ p + dkdidjdiUa,p u*^ p 
+3(9iMQi,pU*j_j, Ua2,pdkdjdiu*^^ p + 



+3 dtdkdlUa^^pU*^^ p Ua2,pdjU*^2,P + 
-1-6 didkUa^,pU*^^ p Ua2.pdjdlU*^^ p], 



and 



Ep hijki{p)r 



(31) 



(32) 
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In the analogy with the Laughlin case we would expect 
that the coefficient of the quartic term is a "square of 
metric" , but the tensor hij^i can not be greatly simplified 
without further assumptions. Assuming Hp = constant, 
we can shift the sum over p to the one over P — p ~ q/2, 



and discuss the product, u* 



this case the quartic coefficient becomes, 



"'ijkl 



fSFS 
Hi 



Qijkl 



separately. In 



(33) 



where Qijui is a tensor. By analyzing u*^ p^^^^Ua^p-q/2 
product and rewriting it as 



^Ip+f "a,P-f = exp{~iq,A^}f{u,u*), 



(34) 



where Ai = —iu'^ pdua^p, the Berry connection, we can 
find out the expression for the gijki tensor in its gauge 
invariant form, 



25. 



ijkl 



{diUa^pdjdkdiu*^ p + didjdkUa,p diu*^ p) 



+ -iAi{ua,pdjdkdiu*^ p - didkdiUa^pU*^^p 

+didkUa,pdjU*^ p - diUa,pdjdiu*^ p) 
+4:A,AjdkUa,pdiul p - 2A,AjAkAi. (35) 

Even the condition on the constancy of the metric, 
didjg^f = 0, does not simplify the form of the gijki ten- 
sor and the coefficient hfj'j^i . If we choose the LLL basis 
of Ref. IM", for which gff = gf/ = B/2, where B is the 
averaged Berry curvature, i.e. 



y B 

Abz 



27rC 



(36) 



where A-qz is the area of the BZ, and gf2 — 92i ^ 



0, 



we expect to recover the "square of metric" form of the 
coefficient, otherwise, for a general FCI, we expect Eq.[33| 
to hold. 



IV. SSF FOR FCI 



in order for the cancelation to occur. 
In the following we will consider that 



B 



9ii = 922 ^ 9>^ and 



9fi = 9if 



0, 



(38) 



i.e. the averaged over BZ FS metric is diagonal with 
diagonal element equal to g. The value of g has the lower 
bound, B/2, as explained in Ref. [T^. The diagonal form 
of the averaged metric will hold in the context of FCI 
states based on the Haldane model Chern insulator, when 
the assumption Up = constant is applied. 

Thus in calculating the SSF for FCI, in the approxima- 
tion we adopted, either we can assume that the system 
in the long-wavelength limit is described by the Laugh- 
lin wave function (or its generalizations with in general 
"incongruent" relationship between flux and particle po- 
sitions or unknown short distance behavior) , and consid- 
ering only the long-wavelength domain apply the plasma 
approach of Section |llj or equivalently we can assume 
that we have a bosonic quantum liquid system with the 
long-range force to which we can apply the Feenberg for- 
mula fin (Eq. [9| in this hmit. 

Either way, assuming the analyticity in m and repeat- 
ing the steps of Eqs. [l6l|T9l we find that the two particle 
part of the SSF of FCI in this limit behaves as 



(9) 



9 n\q[ 



9' n\q\\ 



m B 



-1) 



(39) 



Within the assumptions made (and that we work with 
the FS metric, g[f = g22 = g) we can conclude that the 
form of the projected to a band SSF for a FCI state is 



-FCI 



{<!)= (™^-2) 



-n qtqjqkqi 



9ff9ki^ + 9ijki 



(40) 



A. A mean field picture 

In the following we will consider a possibility that the 
absence of the quadratic term in the expansion of the 
projected SSF in the FQHE also occurs in the context of 
the FCI physics. In a mean field picture we may expect 
that the FCI system in the long-wavelength limit is a 
system with density n — l/27rm/^, fixed by the value of 
= B [15] (i.e. averaged Berry curvature i.e. Chern 
number), and that two body correlations are described 
in the same limit with an effective long-range density- 
density force. 



27rm 



<liqj9.. 



FS ' 



(37) 



In the case g — B /2 we expect to recover the usual QHE 
form. 



B. Discussion 

According to the Ref. |14J fiuctuating Berry curvature 
will produce the quadratic term in the expansion of the fk 
function (Ref. ^) in the SMA for the FCI. The authors 
concluded that necessarily the projected SSF of the FCI 
has to have the leading quadratic term in order to have 
a finite gap in the SMA. 

We applied a mean field approach in calculating two- 
body correlations for the SSF of FCI. It is likely that only 
an exact numerical calculation based on a concrete FCI 
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state may determine whether the quadratic term in the 
expansion of the SSF is present or that the quartic term 
is dominant in determining the physics and energetics of 
FCIs (as in the usual FQHE case). 

If we nevertheless maintain that in a mean field pic- 
ture the formula (40) enters the expression for the gap 
hmction (as a dominator - a norm of the SMA state) in 
the SMA of the FCI state, we can conclude that large 
g B/2 may induce an instability {s^'^^{q) will be- 
come negative, which cannot be true for a positive defi- 
nite quantity) towards a gapless state {s^^^{q) ~ q^). In 
reality we might expect that either large discrepancy be- 
tween g and the lower bound B/2, or strong fluctuations 
of the FCI metric may lead to a gapless state. In order 
to investigate this question we calculated gff , and the 
standard deviation of the FS metric from its averaged 
value, gff ^ in the Brillouin zone for a particular model, 
and we will present this in Section [V] 

To understand better (see also P7H^ ) the absence 
of the quadratic term in the FQHE we will discuss the 
case g — B/2 and _B(Berry curvature) = i? in the FCI 
context. Expanding the expression in Eq. [2T]for the pro- 
jected density to a single band we have for the linear term 
in q: 

Pqllinear = 9^' ^{iAfe(p)7^7p + ^Iq^^p} = l'^^^- (41) 



where : : sign denotes the normal ordering. This defi- 
nition implies the usage of the density operators in the 
SSF calculation that differs by the factor exp{— |gp/2} 
from the usual operators. Nevertheless, this trivial 
difference should not affect the absence of the quadratic 
term, and when applying the expression in Eq. |45[ we see 
that, after a ground state value substraction i.e. normal 
ordering - see Appendix for explanation of this point, the 
quadratic term will not exist if and only if 



(46) 



We know that this is satisfied in the disk and spherical 
geometry of a continuum system [351 130] , and expect to 
hold even in the lattice system, because the generator of 
translation should annihilate the ground state which is a 
homogenous, liquid state. 

Therefore it is the existence of a homogenous ground 
state that is annihilated by the translation generator plus 
the existence of the bosonic algebra of the LLL that en- 
sures the absence of the quadratic term, and appearance 
of the leading quartic term in the FQHE. 

If we believe that the same scenario will happen in 
FCI we may consider the possibility that locally, in the 
BZ, even for varying curvature, we can have the bosonic 
algebra; 



In the first quantization picture the operator Tk is 

N Q 

T, =^{iAfc(K) + ^}, k = x,y (42) 

Using the complex representation we can rewrite the lin- 
ear term in Eq. [41] in the radial gauge as 



The solution must be of the form, 



dp* 



*o = /(fe})exp{-l^|p,p}. 



(44) 



i.e. belong to the LLL, and the operators, Ri = 
(B/4:)p* + d/dp^ and R] = (B/4)pi - 8/ dp*, we rec- 
ognize, corresponding to guiding center coordinates in 
the momentum representation of the QH problem. They 
make the simple bosonic algebra, [ai,al] = 1, if we take 
B = 1 and Ui = \f2Ri and a'^ = \/2RI of the LLL for each 
particle. In this representation SSF can be expressed as 

exp{g*4 - qR^} :: exp{-g*i?] + qR,} :) - 
^(: exp{g*4 - qR,} :) ^(: exp{-g*i?] + qR,} :), 

i j 

(45) 



B + 5B , 

= . Pt 



and i?J 

dpi 



B + 5B d 
4 dp* 



where SB is a weakly dependent function on p. With the 
condition that the FCI state must satisfy. 



E 



i?«*FCl({Pz}) = 0, 



(48) 



and the normal ordering prescription for the bosonic al- 
gebra at each p, the quadratic term will be absent in the 
low-momentum SSF expansion. 

As we will see in a particular example of a bosonic 
FCI state based on the Haldane model at and around 
the point, p = (0, 0) in the BZ, which is the long-distance 
expansion point, in this particular model, the Berry cur- 
vature is zero. Thus the effective form of the ground state 
in this long- wavelength limit (when the particle momenta 
are small) may be of the Jastrow -Laughlin form, i.e. 



'Pj\ 



(49) 



where 7 is a constant, in order to satisfy Eq. [48] in the 
limit for which the Berry curvature and metric matrix 
elements, see Figs. [5| and |6| below, are zero. With the 
assumption that the cancelation of the quadratic terms 
in the small momentum expansion of the projected SSF 
occurs we expect 7 = mB/2g. 
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V. FCIS BASED ON THE HALDANE MODEL 

A. Exact diagonalization results 

The Haldane honeycomb model [T] is the first stud- 
ied example of a Chern insulator. Several numerical ev- 
idences of a robust FCI have been reported for bosons 
with on-site repulsion on such a lattice [9l[3Tl|32]. Weaker 
FCIs have also been observed when bosons are replaced 
by fermions [33] (Note that the on-site interaction is then 
replaced by a nearest neighbor interaction). We will use 
the honeycomb lattice layout of Ref. [2], as shown in 
Fig. [T] The one-body Hamiltonian can be written in 
Bloch form as /i(k) — dol + J^i diCTi using the Pauli ma- 
trices and where 

do ~ 2t2cos0 (cos kx + cos ky + cos{kx + ky)) , 

dx — ii (1 + cos{kx + ky) + cos ky) , (50) 

dy — —t I {8iii{kx + ky) + sin ky) , 

dz = M + 2i2sin0 (sin fc^; + sin fcy — sin(fc3; + /cj,)) . 

M (resp. —M) is the chemical potential added to the 
A (resp. B) sites, ti is the amplitude of the (real) near- 
est neighbor hopping term and ^2 exp(i(/)) the complex 
amplitude if the next nearest neighbor hopping term. 
The two components of the lattice momentum k are de- 
fined as kx = k • ei and ky — k • 62, where ei and 62 
are the lattice vectors. For our numerical calculations, 
we set ti = t2- The Haldane model has two bands. 
If M/ti > 3-\/3 sin((/)), the two bands are trivial. If 
M/ti < 3^/3sin{(j)) then each band carries a non-zero 
Chern number (either C = -|-lorC = — 1). 




FIG. 1: The Haldane model on the honeycomb lattice with 
A (in red) and B (in blue) sublattices. The lattice translation 
vectors are ei and 62. The amplitude of the nearest neigh- 
bor hopping is ti and the next nearest neighbor hopping is 
t2 exp(i(/)) (in the direction of arrows). The sublattice chemi- 
cal potential is set to +M on A sites and —M on B sites. 



We consider N bosons on the Haldane honeycomb 
model with a lattice of unit cells in the ei direction 
and Ny unit cells in the 62 direction. The filling fac- 
tor is thus defined a.s ly — N/{Nx ■ Ny). We add on-site 
Hubbard-type density-density interaction Hiat — J2i ■ 
niUi :, where the sum runs over all the sites. To focus 
on the band topological properties, we use the flat-band 
approximation described in Ref. 4J: We start from the 
original Bloch Hamiltonian /i(k) = X)a=i -^a(k)7-'a(k) 
where -Ea(k) and VaO<-) are the dispersion and the pro- 
jector onto the a-th band, respectively. Then we focus 
on the lowest band and consider the effective one-body 
fiat band Hamiltonian /loff(k) = 7'i(k). From a physical 
perspective, it means that we set the band gap to infinity 
and we make the lowest band completely flat. In this ap- 
proximation, the effective many-body hamiltonian writes 

HcS — 'PlHint'Pl- 



To study the stability of the FCI phase, we focus on 
the energy spectrum. In the FCI regime at filling factor 
z/ = 1/2, the Laughlin-like state on a torus geometry is 
characterized by two almost degenerate low energy states 
separated by a large from higher energy excitations. A 
typical low energy spectrum is shown in the left panel 
of Fig. [2| . The energy splitting between the two lowest 
energy states is called the spread 6. In the case of FQHE, 
the spread should be equal to zero due to the center of 
mass degeneracy. A necessary condition to be able to 
distinguish the two lowest energy states is S to be smaller 
than the gap A (defined as the energy difference between 
the third and the second lowest energy levels). Another 
necessary condition to claim a Laughlin-like is hosted in 
this system is related to the quantum number of the two 
lowest energy states: If they are associated to a Laughlin- 
like state, they should be given by the counting principle 
described in Ref. [T7] . 



We have computed the phase diagram when tuning 
and M at filling ly = 1/2 for two different system sizes: 
TV = 8 on a iV^ = A^j, = 4 lattice (Fig. M and = 10 on 
a Aa; = 5, Ay = 4 lattice (Fig. [s]). We show both the gap 
A and the spread 6 (actually 1 — min((5. A, 1) such that we 
have 1 when the spread is and if (5 > A). When the 
band structure parameters are set to values leading to a 
trivial band, we clearly see that the FCI phase completely 
disappears. In the non-trivial region, both system sizes 
suggest a robust Laughlin-like state around M = and 
(t> = 0.11 - 0.127r. 
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FIG. 2: Left panel: Typical low energy spectrum for the FCI Haldane model as a function of the linearized two-dimensional 
momentum {kx,ky). Here we have set N = 8, Nx = Ny — 4, M = and (p — 0.127r. The spread 5 is the splitting between 
the two lowest energy states, corresponding to the twofold degenerate Laughlin states. The gap A is defined as the energy 
difference between the third lowest energy and the second lowest energy, irrespective of the momentum sector. Middle panel: 
The energy gap A as a function of the two tight-binding model parameters and M/ti, for N = 8 bosons on a A^^; — 4, A^^ — 4. 
The gap is set to zero when the two lowest energy states are not in the expected momentum sectors of the Laughlin state, here 
at {Kx = 0, Ky = 0) for both states. The red line denotes the separation between the Chern insulator phase (upper part) and 
the trivial phase (lower part). Right panel: The corresponding spread 5 (displayed as 1 — mm{S/ A, 1)) as a function of (j> and 
M/ti. A color value of 1 would correspond to a perfectly twofold degenerate groundstate (i.e. 5 = 0). 
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FIG. 3: Left panel: The energy gap A as a function of the 
two tight-binding model parameters (j) and M/t\, for A'' = 10 
bosons on a A^^, = 5, Ny = 4. The gap is set to zero when the 
two lowest energy states are not in the expected momentum 
sectors of the Laughlin state, here at [K^ = 0, Ky — 0) and 
{Kx ~ 0,Ky — 2). The red line denotes the separation be- 
tween the Chern insulator phase (upper part) and the trivial 
phase (lower part). Right panel: The corresponding spread S 
(displayed as 1 — min(5/A, 1)) as a function of (j> and M/ti. 
A color value of 1 would correspond to a perfectly twofold 
degenerate groundstate (i.e. 5 — 0). 



B. The single-particle background 

We now consider the one-body hamiltonian properties. 
In the left panel of Fig. |4j we show the non-universal 



nature of 5 - the averaged over BZ diagonal element of 
the quantum distance (FS) metric, more precisely g^BZ 
(where Abz is the BZ area), is illustrated for the (two 
band) Haldane model with fixed parameters ti ~ t2 ~ 1- 
The standard deviations divided by the averaged values 
g and the standard deviations of the Berry curvatures 
of the Haldane model are shown in the middle and right 
panel of Fig. |4] Note that both the relative deviations 
of the FS metric and those of the Berry curvature are 
minimal around the point (M, (/)) = (0,0.1l7r). This 
is around the same region that we have observed the 
strongest Laughlin-like state in our finite size numerical 
calculations. 



To provide a more detailed insight of this point 
(M, 0) = (0, O.IItt), we provide in Fig. [s] the values of 
the gfj^ tensor (in an orthogonal coordinate system): 

9x1 ,9x2 T 921 1 ^'^d 922 as functions of the Bloch mo- 
mentum in the Brillouin zone. We also give in Fig. |6] the 
values of the Berry curvature at the same point in the 
phase space. 



Discussion plifying procedure of the band flattening the Laughlin- 



In this example of bosonic FCI states based on the 
Haldane model we found that, after the usual but sim- 
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FIG. 4: Left panel: quAbz = qAbz for the Haldane model as a function of (p and M with ti = t2 = 1- Middle panel: The 
corresponding deviation of gn, i.e. Agn/gn. Right panel: The relative deviation of the Berry curvature B with respect to its 
averaged value. 




FIG. 5: The values of the FS metric with respect to its 
mean values at {M,<j)) — (0, O.IItt). The mean values are 
^Bz(5ii) = ABz{g22) = 1.13 and (312) = (321) = in the 
units of Abz|B| — 2n. In the Figure the graphs are ordered as 
the metric matrix elements. All four metric matrix elements 
are zero at k = (0, 0) momentum. 



type ^ 1/2 FCI state with the most significant gap and 
clear two-fold topological degeneracy corresponds to the 
phase-space parameters with the least dispersion of the 
Berry cm'vature and FS metric in the Brillouin zone, near 
FQHE values. Second, along M = line, where we do 
not expect a formation of charge density wave but pos- 
sibly a transition into another liquid state-superfluid, we 
observe the FCI state that is slowly modified but persis- 
tent with increase in angle, the value g of the FS met- 
ric and Berry curvature dispersion. The liquid transition 
should be accompanied with the collapse of a magneto- 
phonon gap and non-analytical be hav ior in the SSF. If we 
approximate the last term in Eq 



40 



with n|q|^(7^/2 and 
then examine the expression for the projected SSF in the 
same equation for the values of g given in the left panel 
of Fig. |4j we see that we are approaching the transition 
around = 7r/2 according to the mean field treatment. 



\B\ 

2 4 6 8 10 



-2|cq 
I 





FIG. 6: The values of the Berry curvature in the BZ at 
(M, (ji) = (0, O.IItt) with respect to its mean value, B = 
— 2n/ABz, denoted by dash line. Plotted values of B are in 
the units of tt/Abz- 



VI. CONCLUSIONS 

Based on the SSF calculations we studied the role of 
band geometry in the context of the FCI physics. We 
identified that the averaged over BZ FS metric plays the 
role of the quantum distance metric in the long-wave 
length domain based on the calculations of the single 
particle part of the projected to the band SSF. We dis- 
cussed the behavior of the complete projected SSF in 
a mean-field framework, and whether and under which 
conditions the scenario of Ref. [32] is possible in the FCI 
context. We illustrated the role of the band geometry 
in the phase diagram of interacting bosons that live on 
the lattice of the Haldane model. The Laughlin 1/ = 1/2 
bosonic FCI state is the most pronounced for the FQHE 
value of the metric, i.e. when the fiuctuations of the FS 
metric and Berry curvature are minimal. The Laughlin- 
like phase persists with the metric increase for a while 
and approaches a transition in an agreement with the 
mean field treatment of the long-distance physics of FCI. 
Further investigations are necessary which may provide 
us also with reasons for occurrence and stability of FCI 
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states. 
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Appendix A: FCI state-single particle contribution 

Here we analyze the expansion in small momentum to 
the quartic order of the expression in Eq. [24] that rep- 
resents the single particle contribution to the projected 
SSF. To fourth order with assumed summations on re- 
peated indices i, j, k 



3! 

and therefore 



didjdkUa,p + o{q'^), 



(Al) 



d-d u* - d d d,u* 

-QidiUa^^p U*^ p -I- qiqjdiUa^,pdjU*^^ p 



HQj Q „, „,* 
2 ' 



2 9idjUa,,p u^^ p - ^qkdidjUa,,p 9fe<,^p 



3! 



-didjdkUa^^p + o{q'^) 



(A2) 



Because u^.pU*^ ^ = 1 and therefore Uc^^pdiu*^ p 
diUa,p w* p = we have to second order 



(P-gP9>|singlc = U + o{q^) . 

To the second order we have 



(A3) 



'^ai,p'^ai,p-qU^^ p_qUa2,p ~ 1 -(- 



[^a,pdi^j'^ct,p ~^ di^j^a,p ^a^p 



-J-2 OiUcK^^p Uq^-^ p '^a2,p^2^ 



'OL2,pl 



Due to 



-^ij[^a,p^i^j^a.,p ~^ ^i^j^a.p ^ot.p 

+2diUa^,pdjU*^^p\ = 



(A4) 



(A5) 



for any symmetric Aij we have 

{P-qPq) Isingic = n - q,qj ^ gfj^{p)np + o{q^) (A6) 

p 

where 

afj^ip) = \[^^Uo,^,p^]K.^,p + djUa,,pdiU*^^ p - 
'^iUai,pUi^^pUa2,pdjU^^p — djUa^ ^pU^^ pUo,2,pdiU^^ p] 

(A7) 

is the Fubini-Study metric. To the third order we find 



^ai ,p'^ai ,p—q'Ua2 ,p—q 



PC, . 



[ua^pdkdidjU*^ p + dkdidjUa^p u^ p 

ididkUai^pU*^^pUa2,pdjU*^^p] + o{q'^) (AS) 

Again differentiating u^^pU*^ p ^ ^ three times we have 

Aijk[ua,pdidjdku*^p + d^djdkUa^p u*^p 
+3didkUai,pdjU*^_^p 

+3diUa,^pdjdkul^p] = (A9) 
and therefore 

''^ai,p^ai,p-qUa2,p-q^oi2,P ~ 1 + QiQj9ij ip) + 

^'^2^'' [dkUa..pdidjU*^ p + d.djUa^pdkul p 

^i^ai ,p^ai ,p^ct2 ,P^^^i^a2 ,p ~^ 
~didkUai.pU^_^pUa2,pdjU^^,p] = 

Mj 



l + q,q^g[S^p) + :i^dk9!,'{p) 



(AlO) 



Appendix B: SSF definition 

By expanding the expression in Eq. [45]we immediately 
see that the quadratic contribution in q will be, 



\q\\RlR,+RM), 



(Bl) 



and it is true that even when the condition in Eq. 46 is 



applied we have a non-zero contribution to the quadratic 
order. The expression in Eq. [45] does not correspond to 
the usual definition of the SSF as in Ref. [211 and only 
after an additional subtraction it reproduces the well- 
known behavior in the classical (continuum) FQHE. 

We will illustrate and explain the source of this discrep- 
ancy in the coordinate representation of the continuum 
FQHE. Fhst, by using the Eqs. [l0]|l4] we see that [22 



%) = .(<?) -n(l-exp{-^}). 



(B2) 



(This is most easily seen in the first quantization picture 
of the many-body problem considering the action of the 
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translation operator in the LLL, exTp{iqd/dz}; we red- 



erived Eq. B2 in Section II A in the second quantization 
as a step towards the discussion of the FCI state.) 

To extract the quadratic contribution to s(q) we con- 
sider the expansion of the unprojected SSF, s{q), and the 
following correlator of unprojected densities, 



h ^{^^{^0. ■ exp{-iq • r^}). 



V 



(B3) 



If we assume the conservation of the angular momen- 
tum and use the complex notation of the LLL for the 
quadratic term, after the subtraction of "self-terms" i.e. 
these generated by X]i(6^P{i*l ' ^i})^ have for the 
quadratic term the following expressions, 



V 4 

1 111'- 



d 



n\q\' 



Zi) 



(B4) 



To get the final expression we used the properties of the 
LLL functions (which are holomorphic up to the Gaus- 
sian factor) and the property of the homogeneity of the 
ground state i.e. that its holomorphic part is annihilated 
by the d/dzi operator (Eq. 46 1 in the momentum rep- 
resentation). 



The substitution of the expression (B4) for s{q) in 



Eq. |B2| would lead to a non-zero, quadratic in q contri- 
bution to s{q). This difference in the value that we have 
for s{q) must stem from a difference in the subtractions: 
the one used in Ref. i22j and the other, when in a static 
correlator the ground state values of the two densities 
are subtracted ("self-terms"), implied by the expression 
of Eq. |45] Namely, in the classical reference the subtrac- 
tion (i.e. a procedure to avoid divergences) is introduced 
for s{q) at any g 7^ as. 



dr(g(r) - 1) exp{-iq • r}. 



(B5) 



because the combination, {g{r) — 1), leads to the ab- 
sence of the divergences for large r. For example, in 
the integer QH case we have (as an exact expression) 
gir) ~ 1 — exp{— |gp/2}, and this leads to the usual, well- 
known behavior s{q) k, n|5p/2 and s{q) ~ 0. In Ref . [22] 
it was shown that s{q) ~ n\q\'^ /2 for any liquid ground 
state of the system that conserves the angular momen- 
tum and particle number. The difference between this 
conclusion and the result in Eq. |B4| stems from different 
subtraction procedures, and can be traced back to two 



different definitions of the SSF, the first in Eq. B5 that 
is related to a static limit of the time ordered density- 
density correlator, and the second, defined as a static 
correlator from which "self-terms" are subtracted (as in 



the expression of Eq. 45 in the momentum representation 
in the projected case). 
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